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Abstract: The sorption-desorption hysteresis observed in many nanoporous solids, at vapor pres- 
sures low enough for the the liquid (capillary) phase of the adsorbate to be absent, has long been 
vaguely attributed to changes in the nanopore structure, but no mathematically consistent expla- 
nation has been presented. The present work takes an analytical approach to account for discrete 
molecular forces in the nanopore fluid and proposes two related mechanisms that can explain the 
hysteresis at low vapor pressure without assuming any change in the nanopore structure. The first 
mechanism, presented in Part I, consists of a series of snap-through instabilities during the filling 
or emptying of non-uniform nanopores or nanoscale asperities. The instabilities are caused by non- 
uniqueness in the misfit disjoining pressures engendered by a difference between the nanopore width 
and an integer multiple of the thickness of a monomolecular adsorption layer. The second mecha- 
nism, presented in Part II, consists of molecular coalescence within a partially filled surface, nanopore 
or nanopore network. This general thermodynamic instability is driven by attractive intermolecular 
forces within the adsorbate and forms the basis to develop a unified theory of both mechanisms. The 
ultimate goals of the theory are to predict the fluid transport in nanoporous solids from microscopic 
first principles, and to determine the pore size distribution and internal surface area from sorption 
tests. 

Introduction 

The sorption isotherm, characterizing the isothermal dependence of the adsorbate mass con- 
tent on the relative vapor pressure at thermodynamic equilibrium, is a basic characteristic of 
adsorbent porous solids. It is important for estimating the internal pore surface of hydrated 
Portland cement paste and other materials. It represents the essential input for solutions of the 
diffusion equation for drying and wetting of concrete, for calculations of the release of methane 
from coal deposits and rock masses, for the analysis of sequestration of carbon dioxide in rock 
formations, etc. Its measurements provide vital information for determining the internal surface 
of nanoporous solids [Ml El EB El EZl SH e.g.]. 

An important feature sorption experiments with water, nitrogen, alcohol, methane, carbon 
dioxide, etc., has been a pronounced hysteresis, observed at both high and low vapor pressures 
and illustrated by two classical experiments in Fig. ^,d) |38, P- 277] and [28j (see also 
[291 Sol El [271 HZl e.g.]). For adsorbates that exist at room temperature in a liquid form, e.g. 
water, the room temperature hysteresis at high vapor pressures near saturation has easily been 
explained by non-uniqueness of the surfaces of capillary menisci of liquid adsorbate in larger 
pores (e.g., the 'ink-bottle' effect [21]). However, a liquid (capillary) water can exist in the 
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Figure 1: (a) Typical desorption and sorption isotherms; (b) BET isotherm; (c)-(d) Desorption and 
sorption isotherms measured on hardened Portland cement paste. 

pores only if the capillary tension under the meniscus (which is given by the Kelvin-Laplace 
equation) does not exceed the tensile strength of liquid water, which is often thought to be 
exhausted at no less than 45% of the saturation pressure, if not much higher. Anyway, at vapor 
pressures less than about 80% of the saturation pressure, the liquid phase represents a small 
fraction of the total evaporable water content of calcium silicate hydrates (C-S-H) [321 Fig. 3] 

However, the hysteresis at low vapor pressures (lower than 80% of saturation in the case of 
C-S-H) has remained a perplexing and unexplained feature for over 60 years. In that case, most 
or all of the adsorbate is held by surface adsorption. The gases and porous solids of interest 
generally form adsorption layers consisting of several monomolecular layers (Fig. [T|d). The 
multi-layer adsorption is described by BET isotherm [22j (Fig. [T^), which is reversible. Sorption 
experiments have generally been interpreted under the (tacit) hypothesis of free adsorption, i.e., 
the adsorption in which the surface of the adsorption layer is exposed to gas. 

In nanoporous solids, though, most of the adsorbate is in the form of hindered adsorption 
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layers, i.e., layers confined in the nanopores (which are sometimes defined as pores < 2nm wide 
[1]). These layers have no surface directly exposed to vapor and communicate with the vapor in 
macropores by diffusion along the layer. It has been well known that a large transverse stress, 
called the disjoining pressure ^] (or solvation pressure [4]), must develop in these layers. 

Development of the theory of hindered adsorption for concrete was stimulated by Powers' 
general ideas on the of creep mechanism [39] . Its mathematical formulation for C-S-H gradually 
emerged in El El E] and was reviewed in a broad context in [9J. But this theory of hindered 
adsorption is also reversible. Thus, although a theory exists, it cannot explain the hysteresis. 

The sorption hysteresis in hardened Portland cement paste, concrete and various solid gels 
[m |33] has for a long time been vaguely attributed to some sort of changes of the nanopore 
structure. In particular, it was proposed that the exit or ingress of water, called the interlayer 
water, from or into the narrowest nanopores would somehow cause large relative changes of 
pore widths ([29j; see also [13], Figs. 13 and 16 in |10], Fig. 1 in [27], or Fig. 9 in [32]), 
picturing pore width changes > 100%. However, if a mathematical model of such a mechanism 
were attempted it would inevitably predict enormous macroscopic deformations, far larger than 
the observed shrinkage caused by drying. Here it will be shown that sorption hysteresis must 
occur even if the nanopore structure does not change. 

The salient feature of capillarity is its non-uniqueness, and the main message of this work 
is that an analogous non-uniqueness also applies to hindered adsorption in nanopores. In the 
capillary range, the non-uniqueness is classically explained by the afore-mentioned "ink-bottle' 
effect, which exists even in two dimensions. In three dimensions, there is much broader range of 
topological and geometrical configurations which provide a much richer and more potent source 
on non-uniqueness of liquid adsorbate content. 

The simplest demonstration is a regular cubic array of identical spherical particles separated 
by a small gap 6 between each pair. At h = 1, either all of the pore space can be filled 
by liquid, or an anticlastic (hyperbolic paraboloid) meniscus surface of zero total curvature 
= rf ^ + = and liquid pressure equal to Ps can exist between each two spheres, with 
ri = — r2 (where ri,r2 = principal curvature radii). This can explain 100% differences among 
equilibrium liquid contents w at h = 1 observed in some experiments. It can even be shown that 
when both S and ri (with r2 = — ri) approach zero in a certain way, then also the liquid content 
(as a continuum) approaches 0. Thus, in theory, an arbitrarily small but nonzero equilibrium 
liquid content at /i = 1 is possible, though extremely unlikely. 

This three-dimensional picture, for example, explains why (as shown in Fig. [l^ by dashed 
lines) the non- uniqueness of sorption isotherm extends to h > 1 (where h = py/ps{T) = relative 
vapor pressure, or relative humidity in the case of water, py = pressure of vapor or gas; Ps{T) = 
saturation vapor pressure). For h > 1, or p^ > ps, the total curvature of the menisci is changed 
from positive to negative, the pores contain overpressurized vapor, and the hysteresis, or non- 
uniqueness continues [121 [13] • This non-uniqueness and hysteresis explains why the slope of the 
isotherm for h > 1 is one, or even two, orders of magnitude higher that one would calculate if 
all the water were liquid for h > 1. (In theory, this nonuniqueness can extend up to the critical 
point of water). In cements these phenomena are complicated by the fact that the chemical 
reactions of hydration withdraw some water from the pores, and create self-desiccation bubbles. 
As a result, one practically never has concrete devoid of any vapor, even for p > p^. 

The consequence of the non-uniqueness is that the sorption isotherm is not a function of 
local thermodynamic variables. Instead, it is a functional of the entire previous history of 
adsorbate content. Here we will show that the same functional character extends to the range 
of hindered adsorption in nanopores, consistent with the extensive experimental data that 
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consistently exhibits sorption hysteresis over the entire range of relative humidities. 

The recent advent of molecular dynamic (MD) simulations is advancing the knowledge of 
nanoporous solids and gels or colloidal systems in a profound way [371 ESI EH ESI EH ESI 1121 
[361 SHI- Particularly exciting have been the new results by RoUand Pellenq and co-workers at 
the Concrete Sustainability Hub in MIT led by Franz- Josef Ulm [HI [TSl US] . These researchers 
used numerical MD simulations to study sorption and desorption in nanopores of coal and 
calcium silicate hydrates. Their MD simulations [HI Fig. 3,4] demonstrated that the filling 
and emptying of pores 1 and 2 nm wide by water molecules exhibits marked hysteresis. 

Especially revealing is the latest paper of Pellenq et al. from MIT [15]. Simulating a 
chain of nanopores, they computed the distributions of disjoining (or transverse) pressure and 
found that it can alternate between negative (compressive) and positive (tensile), depending on 
the difference of pore width from an integer multiple of the natural thickness of an adsorbed 
monomolecular layer (see Figs. 4 and 11 in [IS]. This discrete aspect of disjoining pressure, 
which cannot be captured by continuum thermodynamics, was a crucial finding of Pellenq 
et al. which stimulated the mathematical formulation of snap-through instabilities in Part I 
of this work. Oscillations between positive and negative disjoining pressures have also been 
revealed by density-functional-theory simulations of colloidal fluids or gels in [1] (where the 
excess transverse stress is called the "solvation pressure" rather than the disjoining pressure). 

This work is organized as follows. In Part I, we begin by summarizing the classical theory 
of multilayer adsorption on free surfaces by Brunauer, Emmett and Teller (BET) ^22j, which is 
widely used to fit experimental data, but assumes reversible adsorption without any hysteresis. 
We then develop a general theory of hindered adsorption in nanopores which accounts for crucial 
and previously neglected effects of molecular discreteness as the pore width varies. This leads 
us to the first general mechanism for sorption hysteresis, snap-through instability in nonuniform 
pores, which is the focus of this Part I. 

In Part II [I], we will show that attractive forces between discrete adsorbed molecules can 
also lead to sorption hysteresis by molecular coalescence in arbitrary nanopore geometries, 
including perfectly flat surfaces and pores. This second mechanism for hysteresis is a general 
thermodynamic instability of the homogeneous adsorbate that leads to stable high-density 
and low-density phases below the critical temperature. The mathematical formulation of the 
second part is thus based on non-equilibrium statistical mechanics. Similar models have been 
developed for surface wetting by nanoscale thin films [SH [SSI ISS], starting with Van der Waals 
over a century ago [SZ]- Even more relevant models, accounting for nanoscale confinement, 
have been developed for ion intercalation in solid nanoparticles with applications to Li-ion 
batteries [SH ESI [SSUSO] ■ In that setting, analogous phenomena of hysteresis ^9J (in the battery 
voltage vs. state of charge, in the limit of zero current) and nanoparticle size dependence [H] 
have now been observed in experiments. These connections, which convey the remarkable 
generality of hysteresis in adsorption phenomena, will be developed more in the second part 
in the context of a statistical physics approach. Here, in the first part, we begin to build the 
theory using more familiar models from solid mechanics and continuum thermodynamics. 

Continuum Thermodynamics of Hindered Adsorption in Nanopores 

Free Adsorption: When a multi-molecular adsorption layer on a solid adsorbent surface is 
in contact with the gaseous phase of the adsorbate, the effective thickness a of the layer is well 
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described by the BET equation [22, eq. 28]: 

Sq ii 1 — n 1 — n + ctii 

where T = absolute temperature; = mass of adsorbate per unit surface area; Fi = mass of 
one full molecular monolayer per unit area; G = dimensionless surface coverage; h = relative 
pressure of the vapor in macropores with which the adsorbed water is in thermodynamic equi- 
librium; R = universal gas constant (8314 J kmole"^ °K^^) ; cq = constant depending on the 
entropy of adsorption; AQa = latent heat of adsorption minus latent heat of liquefaction; sq = 
effective thickness of a monomolecular layer of the adsorbate; a = effective thickness of the free 
adsorption layer (in contact with vapor; Fig. [TJj). For the typical value of ct = 54, the BET 
isotherm is plotted in Fig. [TJd, where the number of adsorbed monolayers approaches five at 
the saturation pressure. 

Eq. [T] can be easily inverted: 

h = h{a) = A + VA^ + B (2) 

where A = ^ f 1 - ^ V B = ^ (3) 
2 V a / CT-l ^ ' 

Hindered Adsorption: Consider now a pore with planar rigid adsorbent walls parallel 
to coordinates x and z and a width 2y that is smaller than the combined width 2a of the free 
adsorption layers at the opposite walls given by Eq. ([T]). Then the adsorbate has no surface in 
contact with the vapor and full free adsorption layers are prevented from building up at opposite 
pore walls, i.e., the adsorption is hindered and a transverse pressure, pd, called the disjoining 
pressure [2S], must develop. For water in highly hydrophillic C-S-H, the adsorption layers can 
be up to 5 molecules thick, and so, in pores less than 10 molecules wide {2y < 2.6 nm), hindered 
adsorption with disjoining pressure will develop at high enough h. The adsorbent communicates 
by diffusion of the adsorbate along the pore with the water vapor in an adjacent macropore. 

In a process in which thermodynamic equilibrium is maintained, the chemical potentials /x 
of the vapor and its adsorbate, representing the Gibbs' free energy per unit mass, must remain 
equal. So, under isothermal conditions, 

dfi = p-\dPd + 2dpa)/3 = p-^dp, (4) 

Here p = mass density of the vapor and pa = average mass density of the adsorbate (which 
probably is, in the case of water, somewhere between the mass density pw of liquid water and 
ice). The superior ~ is attached to distinguish the disjoining pressure obtained by continuum 
analysis from that obtained later by discrete molecular considerations {pa = if the nanopore is 
not filled because the transverse pressure due to water vapor is negligible); pa = T^a/y = in-plane 
pressure in the adsorption layer averaged through the thickness of the hindered adsorption layer; 
it has the dimension of N/m^, and (in contrast to stress) is taken positive for compression; tt^ 
= longitudinal spreading 'pressure' in the adsorption half-layer of thickness y (here the term 
'pressure' is a historically rooted misnomer; its dimension is not pressure, N/m^, but force per 
unit length, N/m); tTq is superposed on the solid surface tension gas which is generally larger in 
magnitude, and so the total surface tension, 7 = 7^ — Pa, is actually tensile [71 Fig. 2] (thus the 
decrease of spreading pressure with decreasing h causes an increase of surface tension, which is 
one of the causes of shrinkage). 

Further note that if pd and Pa were equal, the left-hand side would be dp = p~^dpd, which 
is the standard form for a bulk fluid. Also, in contrast to solid mechanics, the left-hand side of 
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Eq. Q cannot be written as tydpd + 2exdpa because strains and cannot be defined (since 
the molecules in adsorption layers migrate and the difference between pd and pa is caused by 
the forces from solid adsorbent wall rather than by strains). 

Consider now that the ideal gas equation PvP^^ = RT/M applies to the vapor (M = 
molecular weight of the adsorbate; e.g., for water M = 18.02 kg/kmole). Upon substitution 
into Eq. Q, we have the differential equation: 

for h < hf p-'dpa = {RT/M) dpjp, (5) 
for h > hf p-\dpd + 2dpa)/3 = {RT/M) dpjp, (6) 

where hf = value of h at which the nanopore of width 2y gets filled, i.e., hf = h{y) based on 
Eq. (|2|. Factors 2 and 3 do not appear ioi h <hf because the free adsorbed layer can expand 
freely in the thickness direction. Integration of Eq. ([6]) under the assumption of constant pa 
yields: 

ioT h<hf: p^ = ^ = Pa^ In h (7) 



for h>hf: Pd + 2{pa - Paf) = ^Pa^^^ 



RT ^ h 
M^^'hf 



where Paf = Pa{hf) = longitudinal pressure when the nanopore just gets filled, i.e., when a = y. 

It is now convenient to introduce the ratio of the increments of in-plane and disjoining 
pressures, 

K = dpa / dpd (9) 

which we will call the disjoining ratio. If the adsorbate were a fluid, k would equal 1. Since 
it is not, K 7^ 1. The role of k is analogous to the Poisson ratio of elastic solids. A rigorous 
calculation of k would require introducing (aside from surface forces) the constitutive equation 
relating pa and pd (this was done in [TT], but led to a complex hypothetical model with too 
many unknown parameters). 

We will consider k as constant, partly for the sake of simplicity, partly because (as clarified 
later) k is determined by inclined forces between the pairs of adsorbate molecules (Fig. [3]d,c); k 
should be constant in multi-molecular layers because the orientation distribution of these forces 
is probably independent of the nanopore width. Note that k would equal only if all these 
intermolecular forces were either in-plane or orthogonal (Fig. [3^, as in a rectangular grid). 

For constant disjoining ratio k, we may substitute pa = upd in Eq. ([8|, and we get 



Pa RT h 
1 + 2k M hf 



For K = 0, this equation coincides with equation 29 in but, in view of Fig. [3| a zero k must 
be an oversimplification. 

According to this continuum model of hindered adsorption, which represents a minor exten- 
sion of [7j, the sorption isotherm of the adsorbate mass as a function of vapor pressure would 
have to be reversible. However, many classical and recent experiments [3H1 EHj SOl ETJ HTJ 
e.g.] as well as recent molecular simulations [HI |T5l [16] show it is not. Two mutually related 
mechanisms that must cause sorption irreversibility in nanopores with fixed rigid walls will be 
presented, one here in Part I, and one in Part II which follows. 
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Fig. 2 

Figure 2: (a) Interatomic pair potential; (b) the corresponding interatomic force and secant stiffness; 
(c) interatomic forces between opposite pore walls visualized by springs. 
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Mechanism I: Snap-Through InstabiUty 

The local transverse (or disjoining) pressure pd can be determined from the transverse stiffness 
Cn, defined as C„ = AF/As where AF transverse resisting force per molecule and As = change 
of spacing (or distance) between the adjacent monomolecular layers in a nanopore containing 
n monomolecular layers of the adsorbate. Since large changes of molecular separation are 
considered, C„ varies with s and should be interpreted as the secant modulus in the force- 
displacement diagram (Fig. |2]d). For this reason, and also because many bond forces are 
inclined ("lateral interactions" [TTllTStfTQ]) 

rather than orthogonal with respect to the adsorption layer (shown by the bars in Fig. [s]), 
Cn is generally not the same as the second derivative d^$/dr^ of interatomic potential nor the 
first derivative dF/ds of force F (Fig. [2^,b). 

To estimate Cn, one could consider various ideahzed arrangements of the adsorbate molecules 
(as depicted two-dimensionally for two different pore widths 2y in Fig. |3| and thus obtain 
analytical expressions for C„ based on the classical mechanics of statically indeterminate elastic 
trusses. However, in view of all the approximations and idealizations it makes no sense to delve 
into these details. 

Diverging Nanopore: Consider now a wedge-shaped nanopore between two diverging 
planar walls of the adsorbent (Fig. |4^), having the width of 2y where y = kx. Here x = 
longitudinal coordinate (Fig. |4]), /c = constant (wedge inclination) and sq = effective spacing 
of adsorbate molecules at no stress. In the third dimension, the width is considered to be 
also Sq. The adsorbate molecules are mobile and at the wide end (or mouth) of the pore they 
communicate with an atmosphere of relative vapor pressure h in the macropores. 

We assume the hindered adsorbed layer to be in thermodynamic equilibrium with the vapor 
in an adjacent macropore. This requires equality of the chemical potentials ft per molecule 
{fi = fi/M, the overbar being used to label a quantity per molecule). At the front of the 
portion of the nanopore filled by adsorbate, henceforth called the 'filling front' (marked by 



circled 2, 3 or 4 in Fig. |4j), Eq. (10) of continuum thermodynamics gives a zero transverse 
pressure, pd = 0, and so Ji = fia = fiv 

However, in the discrete treatment of individual molecules, the chemical potential can be 
altered by transverse tension or compression Apd (Fig. |4]), which can develop at the filling 
front and act across the monomolecular layers unless the nanopore width 2y at the filling front 
happens to be an integer multiple of the unstrained molecular spacing sq. We will call Apd the 
'misfit' (part of) disjoining (or transverse) pressure, by analogy with the misfit strain energy 
for a dislocation core in the Peierls-Nabarro model [20] . 

The misfit pressure, which, at the filling front, represents the total transverse pressure (or 
stress), is determined by the average change As of spacing s between adjacent monomolecular 
layers, which is 

As = 2kx/n-so (n=l,2,3...) (11) 

where n is the number of monomolecular layers across the nanopore width, and sq is the natural 
spacing between the adjacent monomolecular layers in free adsorption, i.e., when the transverse 
stress vanishes (note that for the triangular arrangements in Fig. [2|d,c, Sq is obviously less than 
the natural spacing of unstressed adsorbate molecules, shown as Sq in Fig. |2^). So, the force 
between the molecules of the adjacent layers is F = CAs and the strain energy of the imagined 
springs connecting the molecules is FAs/2 or C(As)^/2 per molecule (if, for simplicity, a loading 
along the secant is considered). 

The hindered adsorbed layer is in a multiaxial stress state, for which the total strain energy 
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Fig. 4 

Figure 4: Filling of a continuously diverging nanopore and disjoining pressures. 
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is the sum of the strain energies of the strain components. Since continuum thermodynamics 
gives zero disjoining (transverse) pressure pd at the filhng front, it suffices to add to C„(As)^/2 
the chemical potential fla per molecule at the filling front due to longitudinal pressure pa only. 



So, in view of Eq. (10), the chemical potential per molecule at the filling front with n 



monomolecular layers is 



f^f,n = -So] +f^a (12) 



n 



where the overbar is a label for the quantities per molecule. Since = at the filling front x*, 
the only source of /i„ is the longitudinal spreading pressure pa in the adsorption layer. 

Let us now check whether at some filling front coordinate x* (Fig. |4]) the diverging nanopore 
is able to contain either n or n + 1 monomolecular layers with the same chemical potential per 
molecule. For n + 1 layers, 

(13) 

Setting fin = fin+i, we may solve for x. This yields the critical coordinate and critical pore 
width for which the molecules in n and n + 1 monomolecular layers have the same chemical 
potential per molecule: 




* _ + yCn+l -So * _ * n/i^i 

^f,n- 1 /r- , 1 /r 9k-' yf,n - '^'^^f,n U^^J 

So the critical relative pore width 2yj / sq at the filling front is a weighted harmonic mean of n 
and 77, + 1 (and a simple harmonic mean if Cn = Cn+i)- 

Equality of the chemical potentials per molecule at the filling front for n and n+1 monomolec- 



ular layers in the same nanopore, which occurs for the pore width given by Eq. (14), implies 
that no energy needs to be supplied and none to be withdrawn when the number of monomolec- 
ular layers is changed between n and n + 1. So the equilibrium content of hindered adsorbate in 
the nanopore for a given chemical potential of vapor is non-unique. Similar to non-uniqueness 
of capillary surfaces, this non-uniqueness underlies the sorption-desorption hysteresis in the 
nanopores. 

Misfit Disjoining Pressure: In view of Eq. (12), its value corresponding to fin for n 
monomolecular layers in the nanopore is 

/ 2kx* \ 

Pd,n = Cniso —j+Pd{Xn) (15) 

where Pd{xn), based on continuum thermodynamics, is non-zero if x„ 7^ Xf,n- In contrast to 
stress, the pressure is considered as positive when compressive. Replacing n with n+1, we 
find that the disjoining pressure makes a jump when the number of monomolecular layers in 
the nanopore changes from to n + 1; 

^Pd,n = Pd,n+1 - Pd,n = 2A;X„ (— - + Sq (C„+i - C„) (16) 

\ n n + 1/ 

(see Fig. |4]). At the filling front, the jump is from transverse tension to compression (Fig. 
[sb). The sudden jumps Apd,n of the misfit pressures from tension to compression diminish with 
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increasing n {n = 1,2,3, ...) as the wedge-shaped nanopore is getting wider; see Figs. ^ and 
4ji. For n > 10, these jumps become insignificant. 

Note that, since the changes As of molecular distance are large, the C values depend on 
F or Apd (Fig. |2|d). So Eq. 16 is actually a nonlinear equation for Apd and its numerical 
solution would require iterations. But here we are aiming at conceptual explanation rather 
than numerical results. 

Misfit Chemical Potentials and Their Effect on Sorption Isotherm: The variation 
of chemical potential at the filling front x/ is shown in Fig. |5ji. Since transverse tension at the 
filling front gives the same chemical potential as transverse compression of equal magnitude, 
the misfit chemical potential, defined as the part of chemical potential due to pa at the filling 
front, varies continuously, provided the pore width varies continuously, too; see Fig. |5]i. This 
is because transverse tension gives the same chemical potential as transverse compression of 
equal magnitude. 

The total chemical potential at the filling front is obtained by adding the chemical potential 
fia{xf) obtained from continuum thermodynamics, which yields the potential variation in Fig. 
5^. Considering the relation of filling front coordinate Xf to the adsorbate mass w shown (in a 
smoothed form) in Fig. [sjo, and the relation h = e^^^^'^'^^f , one can deduce the solid curve in 
Fig. |5^ representing the diagram of equilibrium states of mass content w versus relative vapor 
pressure h in the macropore. 

Why are the segments of the pressure variation in Fig. ^ linear, and why are the segments 
of the chemical potential variation in Fig. |5jl,e,f parabolic? The reason is that the variation 
of nanopore width has been idealized as linear (and that the plots are made for constant C). 
These segments take different shapes for other width variations. 

Sequential Snap-Throughs of Adsorbate Content: In sorption testing and most prac- 
tical problems, the relative vapor pressure h is the variable that is controlled, and the adsorbate 
mass w is the response. Consequently, the states at the reversal points 1, 3 5, 7 of the equilib- 
rium diagram in Fig. |5] for the diverging nanopore are unstable. Likewise the states at points 
1, 3, 5, 7 in Fig. [6jl for the nanopore of step-wise variable width. The loss of stability can 
be evidenced by checking that the molecular potential loses positive definiteness. Fundamental 
though such checks may be, it is simpler and more intuitive to argue in terms of infinitely small 
deviations dh from the equilibrium state. 

Consider, e.g., that, in Fig. [5]F or|6]i, a sufficiently slow gradual increase of h has moved the 
equilibrium state from point 2 to point 3, which is a local maximum of /i as a function of w. For 
a further infinitesimal increase dh there is on the equilibrium diagram no longer any point close 
to point 3. So, borrowing a term from structural mechanics [lOj, we realize that the adsorbate 
mass content w must dynamically 'snap through' at constant h along vertical line 34 to point 
4. After dissipating the energy released along segment 34 (the rate of which depends on the 
lingering times of adsorbed molecules and diffusion along the hindered adsorbed layer [H]), 
thermodynamic equilibrium is recovered at point 4. It is stable because a further infinitesimal 
increment of dh finds, next to point 4, an equilibrium state with adsorbate content incremented 
by dw. 

If h is increased slowly enough further, the equilibrium system will move from point 4 to 
point 5 at which a local maximum of h is reached again and the loss of stability gets repeated, 
since a further increase dh can find equilibrium only after a dynamic snap-through to point 
6. Each snap-through will release some energy which must be damped and dissipated by the 
system. So the local maxima of h at points 1, 3, 5 and 7 are the critical states giving rise to 
the so-called 'snap-through instability' ^lOj. The equilibrium states on curved segments le2. 
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Figure 5: Misfit disjoining pressures and chemical potentials in a continuously diverging nanopore, 
with dynamic snap-throughs of adsorbate content. 
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Figure 6: Step-wise diverging nanopore and misfit disjoining pressures. 
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3c4 and 5a6 are unstable and can never be reached in reality. 

The salient feature is that a different path w{h) is followed when h is decreased. To show 
it, consider point 7 in Fig. [5]F or|6jl as the starting point. During a slow enough decrease h, the 
system will follow the stable states along segment 76a until a local minimum of h is reached 
at point a, which is the stability limit. Indeed, if h is further decremented by dh, there is no 
equilibrium state near point a. So the equilibrium state a is unstable and the system will 'snap 
through' djTiamically at constant h along path ab. At point b stable equilibrium is regained 
after sufficient time. When h is decreased further slowly enough, the equilibrium states move 
through segment Mc until again a local minimum of h is reached and stability is lost at point c. 
Thereafter, the system 'snaps through' along line cd to point d, where equilibrium is regained, 
etc. 

In the diverging pore in Fig. § the snap-through means that when the equilibrium filling 
front reaches the critical points, xi,X2, or X3, it will advance forward a certain distance at 
constant h, as fast as diffusion along the micropore, controlled by the lingering times of the 
adsorbate molecules, will permit. The cross-hatched areas in between the sorption and desorp- 
tion isotherms, such as area SicdS in Fig. |5^ or Fig. [6jl, represent sorption hysteresis. They 
also characterize energy dissipation. 

Sequential Snap-Throughs for Step-Wise Nanopore Width Variation: The dia- 
grams in Fig. |5]i,e,f are valid only for a micropore with continuously diverging rigid planar 
walls (Fig. |5^). This is, of course, an idealization. Because of the atomistic structure of pore 
walls, the pore width in reality varies discontinuously, as exemplified in Fig. [6^. The chance of 
a width exactly equal to an integer multiple of Sq is small. 

Consider that the jumps of nanopore width (Fig. [6^) occur at a;i, X2, xs, and that at 
xi is narrower than sq, at X2 exactly equal to 2so, and at X3 wider than Ssq. Thus the filling 
front in pore segment {xi,X2) is in transverse compression, in segment (0:2, X3) at zero trans- 
verse pressure, and in segment (0:3, X4) in transverse tension; see Fig. [6|d. The corresponding 
strain energies, representing the misfit chemical potential Aft^ per molecule, have a pulse-like 
variation as shown in Fig. Continuum thermodynamics, which ignores the misfits, gives 
a monotonically rising staircase variation of the chemical potential fla{x) (per unit mass) as a 
function of the filling front coordinate Xf, represented by path jiedhgba (Fig. [6]i). Superpos- 
ing on this staircase the misfit chemical potential A/i^ (converted to unit mass), one gets the 
non-monotonic step-wise path of equilibrium states, shown by the bold line 12ed455b6 in Fig. 
ii. 

Taking into account the dependence of the adsorbate mass w in the nanopore on the filling 
front coordinate Xf, one can covert the diagram in Fig. |6]i into the sorption isotherm in Fig. 
[6^, usually plotted as w versus h. The monotonic staircase hkfejicba would represents the 
equilibrium path if the misfit disjoining pressures were ignored. 

When the rise of h, and thus yU/, is controlled, the segments 23 and 56 in Fig. |6^ are unstable 
and unreachable. Indeed, when h or Hf is infinitesimally increased above point 2, there is no 
nearby equilibrium state, and so the system will 'snap through' dynamically to point 3. At that 
point, equilibrium is regained, and h and /i/ can be raised again, slowly enough to maintain 
equilibrium, along path 345. A similar dynamic snap-through is repeated along segment 56, 
after which the stable segment 678 can be followed. Likewise, in the diagram of /i/ versus 
the filling front coordinate Xf (Fig. [6]i), forward snap-throughs at increasing fij (which is a 
monotonically increasing function of h) occur along segments 23 and 56. 

When h 01 IX f is decreased slowly enough from point 8, the stable equilibrium path 8766c is 
followed until stability is lost at point c (Fig. [6fe). Then the system snaps through dynamically 
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from c to d, follows equilibrium path def, and snaps dynamically from / to g. Likewise, in Fig. 
|6]i, backward snap-throughs at decreasing /i/ occur along segments be and ef. 

Obviously, the states on segments c5 and /2 in Fig. |6^, or 2e and 5e in Fig. |6]i, can never 
be reached. They represent unstable equilibrium. The shaded areas glSeg and d46bd represent 
hysteresis, which leads to energy dissipation. 

Snap-Throughs in a System of Nanopores: The diverging nanopore (Fig. |4| [5^ 
and [6^) is not the only pore geometry producing sorption hysteresis. There are infinitely 
many such geometries. In the simple model of discrete monolayers pursued in Part I, the only 
geometry avoiding hysteresis due to sequential snap-throughs is hypothetical — the widths of all 
the nanopores would have to be exactly equal to the integer multiples of the natural spacing sq 
of monomolecular layers in free adsorption, so as to annul the misfit pressures. Below, we will 
show that if molecular coalescence is allowed in the lateral direction, then even these special 
pore geometries will exhibit sorption hysteresis, and so the effect is extremely general. 

An essential feature of nanoporosity is that there are nanopores of many different thicknesses 
2y densely distributed as shown in Fig. [7} At a given vapor pressure, all the nanopores that are 
narrower than a certain width 2y are filled by adsorbed water and the wider ones are empty, 
containing only vapor; see Fig. [7^,c,e. 

As the relative pore pressure h is increased, larger and larger pores fill up. A critical state 
(or a local maximum of h) is reached for a pore width at which the misfit chemical potential 
A/Xrf due to misfit disjoining pressure is for n monomolecular layers equal to or larger than the 
misfit chemical potential for n + 1 layers. After that state, the system loses stability and regains 
it only when all the nanopores up to a certain larger width get filled without increasing h. For 
decreasing h, the stability loss would occur for a different pore width. 

The distribution of nanopore thicknesses 2y may be characterized by a continuous cumu- 
lative frequency distribution function (p{y) that represents the combined volume of all the 
nanopores with thicknesses < 2y. This case, though, is not qualitatively different from the 
diverging nanopore studied previously. For (p{y) oc fcy^, the nanopore system in Fig. [t] becomes 
mathematically equivalent to the linearly diverging nanopore studied before. 

The way the hysteresis in the individual nanopores gets superposed to produce a pronounced 
hysteresis on the macroscale is schematically illustrated in Fig. [8j 

Analogy with Snap-Through Buckling of Flat Arch: There is an instructive analogy 
with the snap-through buckling of elastic arches or shells under controlled load (Fig. [9] [TUl 
p. 231]. If the arch is fiat enough and flexible enough not to fracture, the equilibrium diagram of 
total load p versus midspan deflection u follows the diagram sketched in Fig. |9| The segments 
051 and 432 consist of stable states at which the potential energy is positive definite (i.e., has 
a strict local minimum). But this is not true for the equilibrium states on the segment 12, at 
which the potential energy does not have a strict local minimum. 

Consider that load p is increased from point up to the local maximum at critical point 
1 (Fig. [9]). If load p is increased further by an infinitesimal amount ddp, there is no nearby 
equilibrium state. The arch must follow at constant load the dynamic snap-through path 
14, during which there is accelerated motion, with the load difference from the equilibrium 
curve below being equal to the inertial force, which provides rightward acceleration. The arch 
gains kinetic energy up to point 3, swings over (along a horizontal line), and then vibrates at 
constant load about point 3 until the kinetic energy is dissipated by damping (without damping, 
it would vibrate indefinitely). Then, if the load is increased further, the arch moves through 
stable equilibrium states on the segment 34. 

When the load is decreased, starting at point 4, the arch will follow the stable equilibrium 
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Figure 8: Superposition of hysteretic loops from different nanopores. 
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Figure 10: Energy dissipated by sorption hysteresis on a full h-cyc\e 1 — ^ — )■ 1 (shaded area), 
and dissipation during mid-range reversals (a,b). 



states along segment 432 until a local minimum is reached at point 2. If the load is decreased 
further by an infinitesimal amount dP, there is no equilibrium state near point 2. So the arch 
must snap through dynamically to point 5, the load being again balanced by inertia forces which 
provide leftward acceleration. During this snap-through the arch gains kinetic energy, swings 
over ot the left of point 5 and vibrates about point 5 until the kinetic energy is dissipated by 
damping. Then the load can be decreased further following the stable equilibrium states below 
point 5. 

Note that even though the arch is elastic and the structure-load system is conservative, 
hysteresis is inevitable. During the cycle, the arch dissipates an energy equal to the cross- 
hatched area 51325 in Fig. |4} 

Energy Dissipated by Hysteresis and Material Damage: The Gibbs free energy 
dissipated per unit mass of the nanoporous material is dG = w djj where dfi = ^ dlnh, which 
has in thermodynamic equilibrium the value for the adsorbate species in the vapor and for the 
adsorbed phases. Therefore, the free energy dissipated per unit volume of material due to the 
hysteresis during a complete cycle, e.g., a drying- wetting cycle of hardened cement paste, is 

^G^^l^.lh (17) 

M J h ^ ^ 

Since h is in the denominator, integrability, i.e., the finiteness of AG, requires that lim/i_>o = 
l//i" where < n < 1. Graphically, AG is proportional to the area between the sorption and 
desorption isotherms in the diagram of w/h versus h (Fig. 10). 
The energy could be dissipated in two ways: 
1) By internal friction in the adsorbed fluid during the dynamic snap-throughs (or the 
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molecular coalescence phenomena discussed in Part II [T]), or 

2) by fracturing or plastic damage to the nanopore surfaces. 
However, the latter seems unlikely since it could be associated with every disjoining pressure 
change and not particularly with the snap-through. The existence of the former is undeniable, 
and the point here is to show that the hysteresis is perfectly explicable without postulating any 
damage to the nanopore surface. 

Anyway, the degree of material damage due to a drying-wetting cycle, if any, could be 
checked by measuring the strength or the fracture energy, or both, of the material before and 
after the cycle. This would have to be done slowly enough on thin enough specimens having 
drying half-times less than 1 hour (< 1 mm thick for cement paste), in which the relative 
humidity h in the capillary pores can be changed without creating a significant gradient of h 
across the specimen wall (in thicker samples, most of the material damage is done by non- 
uniform shrinkage stresses engendered by non-uniformity of h across the wall thickness ^5]). 
Shrinkage and creep experiments on such specimens have been performed at Northwestern [H] , 
but no cycles were performed and strength changes were not checked. It could also be checked 
whether the snap-throughs might be associated with the acceleration of concrete creep due to 
simultaneous drying, called the drying creep (or Pickett effect). 

Sorption Potential: Note that, based on the derivation of Eq. (17), it further follows 
that 

RT w dG 

In other words, the Gibbs's free energy per unit mass of adsorbate as a function of /i is a 
potential for the adsorbate content parameter (3 during a one-way change of h. 



Conclusions of Part I 

We can summarize the findings of the first part as follows: 

1. One mechanism that must be causing sorption hysteresis at low vapor pressure is a series 
of snap-through instabihties causing path-dependent non-uniqueness of adsorbate content 
and dynamic jumps of water content of nanopores at constant vapor pressure. 

2. The snap-through instabilities are a consequence of the discreteness of the adsorbate, 
which leads to non-uniqueness of mass content and to misfit disjoining (transverse) pres- 
sures due to a difference between the pore width and an integer multiple of the thickness 
of a transversely unstressed mono molecular layer of the adsorbate. 

3. The hysteresis is explained by the fact that the snap-through instabilities for sorption and 
desorption follow different paths. 

4. The snap-through instabilities are analogous to snap-through buckling of arches and shells, 
long known in structural mechanics. They cause hysteresis and energy dissipation even 
when the arch or shell is perfectly elastic. 

If a quantitative version of this theory were developed, it might be possible to infer from the 
hysteresis the surface area and the size distribution of the nanopores filled by hindered adsor- 
bate. Our preliminary analysis of snap-through instabilities suggests that the key to making this 
connection is to account for inclined forces, or "lateral interactions" , in the statistical thermo- 
dynamics of hindered adsorption. In the Part II, we will show that attractive lateral interactions 
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generally lead to sorption hysteresis in any pore geometry due to molecular coalescence of the 
adsorbate. 

Acknowledgment: The research was funded partly by the U.S. National Science Foundation under Grant 
CMS-0556323 to Northwestern University (ZPB) and Grant DMS-0948071 to MIT (MZB) and partly by 
the U.S. Department of Transportation under Grant 27323 provided through the Infrastructure Technology 
Institute of Northwestern University (ZPB). Thanks are due to Franz-Josef Ulm and Rolland J.-M. Pellenq 
of MIT for stimulating discussions of disjoining pressure based on MD simulations, and to Laurent Brochard 
and Hamlin M. Jennings for further valuable discourse. 

References 

[1] Bazant, M. Z. and Bazant, Z. P. "Theory of sorption hysteresis in nanoporous solids: II. Molecular 
coalescence" . preprint. 

[2] Adolphs, J., and Setzer, M.J. (1996). "A model to describe adsorption isotherms." J. of Colloid 

and Interface Science 180, 70-76. 
[3] Adolphs, J., Setzer, M.J. and Heine, P. (2002). "Changes in pore structure and mercury contact 

angle of hardened cement paste depending on relative humidity." Materials and Structures 35, 

477-486. 

[4] Balbuena, P.B., Berry, D., and Gubbins, K.E. (1993). "Solvation pressures for simple fluids in 

micropores." J. Phys. Chem. 97, 937-943. 

[5] Bazant, Z.P. (1970). "Constitutive equation for concrete creep and shrinkage based on thermo- 
dynamics of multi-phase systems." Materials and Structures, 3, 3-36 (reprinted in Fifty Years of 
Evolution, of Science and, Technology of Building Materials and Structures, Ed. by F.H. Wittmann, 
RILEM, Aedificatio Publishers, Freiburg, Germany 1997, 377-410). 

[6] Bazant, Z.P. (1970). "Delayed thermal dilatations of cement paste and concrete due to mass 
transport." Nuclear Engineering & Design, 24, 308-318. 

[7] Bazant, Z.P. (1972). "Thermodynamics of interacting continua with surfaces and creep analysis 
of concrete structures." Nuclear Engineering and Design, 20, 477-505. 

[8] Bazant, Z.P. (1972). "Thermodynamics of hindered adsorption with application to cement paste 
and concrete." Cement and Concrete Research, 2, 1-16. 

[9] Bazant, Z.P. (1975). "Theory of creep and shrinkage in concrete structures: A precis of recent 
developments". Mechanics Today, ed. by S. Nemat-Nasser (Am. Acad. Mech.), Pergamon Press 
1975, Vol. 2, pp. 1-93. 

[10] Bazant, Z.P., and Cedolin, L. (1991). Stability of Structures: Elastic, Inelastic, Fracture and 
Damage Theories, Oxford University Press, New York; 2nd. ed. Dover Publ. 2003; 3rd ed. World 
Scientific Publishing, Singapore-New Jersey-London 2010. 

[11] Bazant, Z.P., and Moschovidis, Z. (1973). "Surface diffusion theory for the drying creep effect in 
Portland cement paste and concrete" J. Am. Ceramic Soc, 56, 235-241. 

[12] Bazant, Z.P., and Thonguthai, W. (1978). "Pore pressure and drying of concrete at high temper- 
ature." Proc. ASCE, J. of the Engrg. Mech. Div., 104, 1058-1080. 

[13] Bazant, Z.P., and Kaplan, M.F. (1996). Concrete at High Temperatures: Material Properties 
and Mathematical Models, Longman ( Addison- Wesley) , London; 2nd ed. Pearson Education, 
Edinburgh, 2002. 

[14] Bonnaud, P.A., Coasne, B., and Pellenq, R.J.-M. (2010). "Molecular simulation of water confined 

in nanoporous silica." J. Phys.: Condens. Matter, 284110 (15pp). 
[15] Brochard, L., Vandamme, M. and Pellenq, R.J.-M. (2011). Competitive adsorption of carbon 

dioxide and methane in coal." Privately communicated manuscript, submitted to PNAS. 
[16] Brochard, L., Vandamme, M. and Pellenq, R.J.-M. (2011). Poromechanics of nanoporous media." 

Privately communicated manuscript; submitted to J. of the Mech. and Phys. of Solids. 



20 



[17] Nikitas, P. (1996). "A simple statistical mechanical approach for studying multilayer adsorption: 
Extensions of the BET adsorption isotherm." J. Phys. Chem. 100, 15247-15254. 

[18] Cerofolini, G.F., Meda, L. (1998). "A theory of multilayer adsorption on rough surfaces in terms 
of clustering and melting BET piles." Surface Science 416, 402-432. 

[19] Cerofolini, G.F., Meda, L. (1998). "Clustering and melting in multilayer equilibrium adsorption." 
J. of Colloid and Interface Science 202, 104-123. 

[20] Hirth, J. P. and Lothe, J., Theory of Dislocations, (Kreiger, 1992). 

[21] Brunauer, S. (1943). The adsorption of gases and vapors. Princeton University Press, Princeton, 
NJ, 1943 (p. 398). 

[22] Brunauer, S., Emmett, P.T., and Teller, E. (1838). "Adsorption of gases in multi-molecular 

layers." J. Amer. Chemical Soc. 60, 309-319. 
[23] B. Coasne, A. Galarneau, F. Di Renzo, and R. J.-M. Pellenq (2009). "Intrusion and retraction of 

fluids in nanopores: effect of morphological heterogeneity", J. Phys. Chem. C 113, 1953-1962. 
[24] B. Coasne, A. Galarneau, F. Di Renzo, and R. J.-M. Pellenq (2008). "Molecular simulation of 

adsorption and intrusion in nanopores." Adsorption (2008) 14: 215-221. 
[25] B. Coasne, F. Di Renzo, A. Galarneau and R.J.-M. Pellenq (2008). "Adsorption of simple fluid on 

silica surface and nanopore: effect of surface chemistry and pore shape." Langmuir 24, 7285-7293 
[26] Derjaguin, B.V. (1940). "On the repulsive forces between charged colloid particles and the theory 

of slow caogulation and stability of lyophole sols. Trans, of the Faraday Society 36, 203, 730. 
[27] Espinosa, R.M., and Franke, L. (2006). "Influence of the age and drying process on pore structure 

and sorption isotherms of hardened cement paste." Cement and Concrete Research 36, 1969-1984 

(Figs. 2, 6, 16). 

[28] Feldman, R.F., and Sereda, P.J. (1964). "Sorption of water on compacts of bottle hydrated 
cement. I: The sorption and length-change isotherms." J. Appl. Chem. 14, p. 87. 

[29] Feldman, R.F., and Scrcda, P.J. (1968). "A model for hydrated Portland cement paste as deduced 
from sorption-length change and mechanical properties." Materials and Structures 1 (6) (1968) 
509-520. 

[30] Frenkel, D., and Smit, B. (1996). Understanding Molecular Simulations. Academic Press, Sand 
Diego-New York. 

[31] Jennings, H.M. (2000). "A model for the microstructure of calcium silicate hydrate in cement 

paste." Cement and Concrete Research 30, 101-116. 
[32] Jennings, H.M. (2010). "Pores and viscoelastic properties of cement paste." submitted to an 

Elsevier Science journal. 

[33] Jennings, H.M., BuUard, J.W., Thomas, J.J., Andrade, J.E., Chen, J.J., and Scherer, G.W. 

(2008). "Characterization and modeling of pores and surfaces in cement paste: Correlations to 

processing and properties." J. of the Advanced Concrete Technology 6 (1), 5-29. 
[34] Jonson, B., Nonat, A., Labbez, C, Cabane, B., and Wennerstrom, H. (2005). "Controlling the 

cohesion of cement paste." Langmuir 21, 9211-9221. 
[35] Jonson, B., Wennerstrom, H. Nonat, A., and Cabane, B. (2004). "Onset of cohesion in cement 

paste." Langmuir 20, 6702-6709. 
[36] Malani, A., Ayappa, K.G., and Murad, S. (2009). "Influence of hydrophillic surface specificity on 

the structural properties of confined water." J. Phys. Chem. 113, 13825-13839. 
[37] Pellenq, R. J.-M., Kushima, A., Shashavari, R., Van Vliet, K.J., Buehler, M.J., Yip, S. and Ulm, 

F.-J. (2010). "A realistic molecular model of cement hydrates." Proc. Nat. Academy of Sciences 

106 (38), 16102-17107. 

[38] Powers, T.C., and T.L. Brownyard (1946), "Studies of the Physical Properties of Hardened 
Portland Cement Paste. Part 2. Studies of Water Fixation". J. of the Am. Concrete Institute 
18(3), 249-336. 

[39] Powers, T.C. (1966). "Some observations on the interpretation of creep data." Bulletin RILEM 
(Paris), p. 381. 



21 



Rarick, R.L., Bhatty, J.W. and Jennings, H.M. (1995). "Surface area measurement using gas 
sorption: Application to cement paste." Material Science of Concrete IV, editied by J. Skalny 
and S. Mindess, Am. Ceramic Soc, Chapter 1, pp. 1-41. 

Scherer, G.W. (1999). "Structure and properties of gels." Cement and Concrete Research 29 
(1999), 1149-1157. 

Smith, D.E., Wang, Y., Chaturvedi, A., Whitley, H.D. (2006). "Molecular simulations of the 
pressure, temperature, and chemical potential dependencies of clay swelling." J. Phys. Chem. B 
110, 20046-20054. 

Thomas, J.J., Allen, A.J., and Jennings, H.M. (2008). "Structural changes to the calcium hydrate 
gel phase of hydrated cement with age, drying and resaturation." J. of the Am. Ceramic Soc. 91 
(10), 3362-3369. 

Bazant, Z.P., Asghari, A. A., and Schmidt, J. (1976). "Experimental study of creep of hardened 
cement paste at variable water content." Materials and Structures (RILEM, Paris), 9, 279-190. 
Bazant, Z.P., and Raftshol, W. J. (1982). "Effect of cracking in drying and shrinkage specimens." 
Cement and Concrete Research, 12, 209-226; Disc. 797-798. 

M. Vandamme, L. Brochard, B. Lecampion and O. Coussy (2010) "Adsorption and strain: The 
C02-induced swelling of coal," J. of the Mechanics and Physics of Solids 58, 1489-1505. 
Baroghel-Bouny, V. (2007) "Water vapour sorption experiments on hardened cementitious ma- 
terials Part I: Essential tool for analysis of hygral behaviour and its relation to pore structure" 
Cement and Concrete Research 37, 414437. 

Wagemaker, M., Singh, D. P., Borghols, W. J. H., Lafont, U., Haverkate, L., Peterson, V. 
K. and Mulder, F. M. (2011) "Dynamic Solubility Limits in Nanosized Olivine LiFeP04", 
J. Am. Chem. Soc, 133, 1022210228. 

Dreyer, W., Jamnik, J., Guhlke, C, Huth, R., Moskon, J., and Gaberscek, M. (2010). "The 
thermodynamic origin of hysteresis in insertion batteries." Nature Materials 9, 448-451. 
Dreyer, D., Guhlke, C., and Huth, R. (2011) "The behavior of a many-particle electrode 
in a lithium-ion battery", Physica D 240, 10081019 

G. K. Singh, G. Ceder and M. Z. Bazant (2008) "Intercalation dynamics in rechargeable 
battery materials: General theory and phase transformation waves in LiFeP04," Elec- 
trochimica Acta 53, 7599. 

D. Burch and M. Z. Bazant (2009) "Size dependent spinodal and miscibility gaps in 
nanoparticles," Nano Letters 9, 3795. 

Bai, P., Cogswell, D. A., and Bazant, M. Z. (2011) "Suppression of phase separation 
in LiFeP04 nanoparticles during battery discharge", preprint, arXiv: 1108. 2326 vl [cond- 
mat.mtrl-sci] 

de Gennes, P. G. (1985) "Wetting: statics and dynamics", Rev. Mod. Phys. 57, 827-863. 
Gouin H. and Gavrilyuk, S. (2008) "Dynamics of liquid nanofilms", Int. J. Eng. Sci. 46, 
1195-1202. 

Gouin H. (2009) "Liquid nanofilms: A mechanical model for the disjoining pressure ", Int. 
J. Eng. Sci. 47, 691-699. 

van der Waals, J. D. (1893). The thermodynamic theory of capillarity under the hypoth- 
esis of a continuous variation of density. (Translation by J. S. Rowlinson (1979)) Journal 
of Statistical Physics, 20, 197 (Original version: Zeitscherift fuer Physikalische Chemie, 
Stoechiometrie und Verwandtschaftslehre, 13, 657). 



22 



